Disorder and transport in cuprates: weak localization and magnetic contributions 
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We report resistivity measurements in underdoped YBa2Cu306.6 and overdoped Tl2Ba2Cu06+a; 
single crystals in which the concentration of defects in the Cu02 planes is controlled by electron 
irradiation. Low T upturns of the resistivity are observed in both cases for large defect content. In 
the Tl compound the decrease of conductivity scales as expected from weak localization theory. On 
the contrary in YBa2Cu306.6 the much larger low T contribution to the resistivity is proportional 
to the defect content and might then be associated to a Kondo like spin flip scattering term. This 
would be consistent with the results on the magnetic properties induced by spinless defects. 
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The origin of the strange metallic properties of the 
cuprates, among which the transport properties, is still 
the key feature to be resolved to progress in the under- 
standing of the physics of these 2D correlated electron 
systems. As the pseudogap is evidenced both in mag- 
netic and transport properties, some interplay between 
spin and charge degrees of freedom is required. This in- 
terplay can be studied as well from the influence of de- 
fects and impurities on the transport properties. Point 
defects such as Zn impurities induce a strong depression 
of Tc as expected for d-wave superconductivity. This is 
often considered to result from pure potential (unitary?) 
scattering However it has been shown that such 

spinless defects induce a local moment in their vicinity 
|,|. One can therefore wonder whether these magnetic 
degrees of freedom play a role in the charge scattering, 
and induce for instance a (Kondo like) spin flip contribu- 
tion to the resistivity. 

In presence of Zn impurities, low T upturns of the 
resistivity with a InT variation have been observed in 
La2Sr2Cu04 (LSCO) and YBazCuaOr-i (YBCO) 
However these observations have required large Zn con- 
tents or low hole doping in order to suppress supercon- 
ductivity. Data were then limited to large values of the 
2D resistivity which approach or exceed ~ bkfl/D for 
which deviations from weak localization predictions start 
to be observed in 2D systems . These experimental 
constraints have up to now prevented any thorough study 
either of the impurity contributions to the scattering, or 
of the applicability of 2D weak localization theory |^2| . 

Recent accurate NMR measurements in the case of Li 
impurities, which behave similarly to Zn, have evidenced 
a Kondo like variation of the induced local moment sus- 
ceptibility with a Kondo temperature Tk which increases 
with hole doping Q. One might expect that magnetic 
contributions to the transport properties should then dif- 
fer markedly in underdoped and overdoped systems while 
localization effects should always be present. 

In this work, we address both questions for the first 
time by comparing underdoped and overdoped cuprates. 
We have used electron irradiation to introduce in-plane 



point defects which behave similarly to Zn substitutions 
in YBCO [p"3p|. As measurements could be performed 
on a single sample, the accuracy on the defect contri- 
bution to the resistivity was high enough to study thor- 
oughly its variation with defect content and tempera- 
ture, even when Tc is still high. In the heavily overdoped 
Tl2Ba2Cu06+a; (T12201), we find here that can be sup- 
pressed while the low T resistance is still below the strong 
localization limit. This enables us to demonstrate that 
weak localization theory applies and that defect scatter- 
ing is purely elastic in this system. We shall show that 
this approach fails in underdoped YBCOg.e- The large 
magnitude of the low T upturns of the resistivity and 
its variation with defect content rather indicates the oc- 
curence of a Kondo like inelastic spin-flip scattering. The 
difference between these two systems will be found con- 
sistent with the observed doping dependence of Tk ■ 

We have taken data on the underdoped YBCOg.e sin- 
gle crystal, which has already been used in ref Q| . It dis- 
played the usual pseudo gap resistivity temperature de- 
pendence with an initial Tc of 60K. The highly overdoped 
T12201 sample prepared similarly to those of ref with 
an initial Tc of 30K, exhibited a power law dependence 
of p {T) characteristic of the overdoped state . The 
samples have been irradiated with 2.5 MeV electrons'^*) 
in liquid hydrogen to allow high electron flux without 
sample heating. After a short time defect annealing (w 
10 minutes) either at ra=100K, 150K or 300K, resistiv- 
ity measurements were performed between 20K and Ta- 
in such conditions the resistivity versus T curves plot- 
ted in Fig.l were found absolutely reversible, indicating 
the absence of further evolution of the defect content. 
The high fluences required to reach large enough resis- 
tivities led us to perform several independent irradiation 
runs. This forced us to remove the probe from the irra- 
diation facility and to let the samples stay at room T 
for long time (typically a few months) after each run, 
which reduced subsequently the concentration of defects 
Q. The resistivity curves after annealing at 300K are 
quite identical to those observed in situ for lower irradia- 
tion fluences. This indicates that the measured resistivity 
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is only sensitive to the concentration of in-plane defects. 
After the last run the samples were removed from the 
irradiation probe and ex-situ measurements were then 
performed down to lower T (1.7K) in a different set-up 
in which magnetoresistance data could be taken. 

As can be seen in Fig.l the high T parts of the resistiv- 
ity curves are nearly parallel to each other in YBCOe.e, 
i.e. Mathiessen's rule extends to large irradiation flu- 
ences Q, which implies also that the hole doping of the 
Cu02 planes is not significantly modified. Conversely the 
variation with T of the resistivity of the T12201 sample in- 
creases with defect content showing that in this case the 
hole doping increases sHghtly [|l5|. In both compounds 
low T upturns of the resistivity are observed and increase 
with defect content. These upturns appear at higher T 
for YBCOg.e, and correspond to larger contributions to 
the 2D in-plane resistance i?2D than for T12201 

To see this effect more quantitatively, we have re- 
ported in Fig. 2 6R2D, the low T increase of i?2_D obtained 
in these two compounds for comparable residual resis- 
tances i?2_D — 4.8fcr2/n. In order to evaluate 5R2d 
the resistivity of the non irradiated compound had to 
be extrapolated below Tc- For pure YBCOg.e, two dif- 
ferent empirical forms R2d{T) = a -\- bT + cT^ and 
i?2£) (T) = a + bT^ have been used to fit the data between 
80K and 150K. This lower T limit has been taken to avoid 
the T range for which superconducting fluctuations con- 
tribute to the resistivity. Although their extrapolations 
differ significantly at low T, the magnitude and T depen- 
dence obtained for SR2D are quite identical in the limited 
experimental range above Tc. For Tl 2201, superconduc- 
tivity has been suppressed for this fluence and the low T 
extrapolation might be more critical than for YBCOe.e- 
To take into account the slight change of hole doping 
we have then fitted the high T part of each resistivity 
curve with R^oiT) = a + bTP. To test the validity of 
this procedure, two extreme cases have been considered. 
The fits were performed either between 70K and 150K, 
or between 150K and 300K. Although SR2D is naturally 
found to depart from zero at a slightly higher T in the 
latter case, the overall variation and particularly the In T 
dependence found at low T is not sensitive to the func- 
tion chosen for i?2£)(T). It is clear in Fig. 2 that, whatever 
the fitting procedure, 6R2D is larger and has a steeper 
T variation in YBCOg.e than in the Tl compound. 

We can at this stage compare these data with the ex- 
pectations of the weak 2D localization theory, which can 
apply as long as the elastic mean free path I verifies 
kpl >> 1 (where kp is the Fermi wave vector), that 
is for i?2D << 26fcr2/n. The weak localization contribu- 
tions are expected to become sizable at temperatures for 
which the elastic scattering rate 1 /t becomes much larger 
than the inelastic scattering rate l/r^. Assuming that the 
in-plane defects only contribute to the elastic scattering, 
the ratio r^/r can be estimated from the data, if the 
inelastic rate is associated to the T dependent contribu- 
tion to R2D of the starting pure material. For instance. 



for the residual resistivity i?2D — 4.8fcJ7/n of Fig. 2, r^/r 
w 10 - 12 and « 7 respectively for T12201 and YBCOe.e 
at lOOK ( at 50K the deduced value would be n/r «40 for 
the two systems). Therefore, within the above assump- 
tions, weak localization conditions are fulfilled even at 
such high T. Let us recall that the theoretical correction 
to the 2D conductance has a universal behaviour at low 
T given by 

= ^p\n(T/To) (1) 

where p is the exponent of the T dependence of the in- 
elastic scattering rate l/r^ cc , and a = 1 in the ab- 
sence of electronic interactions, in good agreement with 
experimental observations |l7|-pOt . 

We have plotted in Fig. 3a the value of Act = 
{R2d)^^ — {R2d)~^ ^'^^ the Tl compound. For the in-situ 
measurements the data are limited to 20K, while they 
could be pursued down to 3K for the final state of the 
sample. In this case a good fit with a InT dependence 
extends over a decade. One can notice in Fig. 3a that the 
curves for the other defect contents are parallel to each 
other at 20K, with a slope identical to that found for the 
lower T run, within experimental accuracy. Finally at a 
given T the magnitude of Act increases with defect con- 
centration, as the upturn in the resistivity curve occurs at 
higher T as can be seen in Fig. la. These qualitative ob- 
servations are more striking in Fig. 3b where we have per- 
formed the scaling of these curves in ln(r/ro). Here Tq 
is found to increase progressively from 70K to 120K with 
defect content. The actual slope obtained experimentally 
below 30K yields the numerical value ap = 1.3, when 
compared to (|l|). The value of p can be independently 
estimated from the T dependence of i?2r> (^) and varies 
from 1.7 to 2.3. This would correspond to a ~ 0.7. These 
results provide a strong experimental evidence that 2D 
weak localization theory explains the T12201 data. 

Transverse magnetoresistance measurements per- 
formed up to ST after the last irradiation provide an 
independent evidence which confirms this analysis. It 
indeed exhibits a negative contribution which decreases 
with increasing T as expected for weak localization. The 
field variation was accurate enough at low T to allow an 
analysis along the lines of ref. |lj] , which allowed to de- 
duce a ~ 0.5 and r^/r ~ 60 at 4.2 K. While the value of 
a is in good agreement with the T dependence of R2D, 
the value of /r is smaller than that which would be ob- 
tained from the extrapolation of the T dependent part of 
the resistivity of the pure sample. Such a trend is quite 
common, as is usually found to fiatten at low T, at val- 
ues smaller, even in order of magnitude, than expected 
from the pure systems estimates ]l7| , p^pll . Further ex- 
periments are needed to determine whether interaction 
effects contribute to Act . 

Let us consider now the data for the underdoped com- 
pound YBCOg.e, for which we have seen that r^/r is 
weaker or similar, whereas SR2D is larger than for the 
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overdoped compound. In this case superconductivity pro- 
hibits to investigate the low T range, and weak locahza- 
tion contributions are not expected to be fully developed. 
However, as can be anticipated from Fig. 2, a plot simi- 
lar to that of Fig. 3 gives much larger values of Act than 
those found for T12201. The slope found at the lowest T 
for which T/Tq is only 0.5 would extrapolate at low T to 
a value ap > 5, which is absolutely not sound physically. 
We could obtain a similar conclusion for electron irradi- 
ated YBCO7 for which p = 1 for the well known T linear 
dependence of p(T) . Therefore all these data cannot 
be explained by weak 2D localization theory, as was also 
concluded for Zn substitution in YBCOg.y 0. 

Alternatively, we can wonder whether spin flip scat- 
tering effects can account for these large low T contri- 
butions. While for weak localization 5R2d scales with 
from (|l|), a Kondo like contribution to R2D is ex- 
pected to scale linearly with the defect content, that is 
with i?^, the increase of the residual resistivity. Indeed 
one can see in Fig. 4 that 6R2D/Rd follow the same T 
dependence for a large range of Rd values, if one omits 
to consider the dowturn associated with the proximity of 
the superconducting transition. Upward deviations with 
respect to this common behaviour begin to occur at low 
T ior Rd > 3.8kfl, that is R2D ^ 5fcf2. This is even more 
prominent for the largest defect content in Fig. lb which 
varies faster than InT, and might therefore correspond 
to the onset of strong localization. 

It is interesting at this stage to compare these results 
with that obtained for underdoped Zn substituted sam- 
ples 0. In that case refined analyses, as done here, are 
impossible in the high T range as the substituted sample 
has to be compared to a pure sample which is physi- 
cally distinct. Systematic errors due to differences in the 
sample geometries and hole contents can then be hardly 
avoided. However superconductivity could be nearly sup- 
pressed with a 18 T field in YBCOe.r with 2.7%Zn @] 
which allowed to evidence a InT contribution to the re- 
sistivity (not to the conductivity) below 30K. This would 
be compatible with a Kondo contribution with a low Tk 
value. A rough estimate for 6R2D/Rd yields a value for 
this InT contribution about a factor two smaller than 
that obtained here. This order of magnitude agreement 
is quite satisfactory as we are comparing here the 
influences of different microscopic defects. 

In conclusion, we have evidenced that the low T resis- 
tivity behaves quite differently in overdoped and under- 
doped systems. In the former case we have shown the rel- 
evance of weak localization theory, as Act has a In T vari- 
ation independent of defect content and in rough quanti- 
tative agreement with the exponent of R2d{T). A nega- 
tive transverse magnetoresistance has been evidenced at 
low T, but further experiments are needed to determine 
the eventual contribution of interactions to the In T term 
and the low T variation of . 

We have shown that the large SR2D measured in the 
underdoped case can be interpreted in terms of a Kondo 
scattering on the local moment induced by the "spin- 



less" defects. This implies per se that the scattering by 
such defects should saturate at low T [|5| and should 
only become unitary below Tk- Therefore, pure poten- 
tial scattering by spinless impurities is only expected for 
T << Tk- As large Tk values, exceeding 200K, have 
been evidenced by NMR only for the overdoped case [|| , 
pure elastic scattering is expected at low T for these sys- 
tems. This is then quite coherent with the occurence of 
weak localization contributions in the T12201 compound, 
if the NMR observations for Li impurities do apply as 
well for irradiation defects. 

It is worth pointing here that -InT contributions to 
R2D were also observed for the "pure" underdoped LSCO 
|6| and Bi(SrLa)CO 2201 They have been asso- 

ciated with a "peculiar" charge localization leading to 
a metal insulator transition (MIT) with decreasing hole 
content. It is usually assumed that some disorder is re- 
sponsible for the anomalously low value of the optimum 
Tc in these compounds and that this disorder might be 
at the origin of the -InT term. Our data lead us to 
suggest that this unknown disorder might result primar- 
ily in a scattering on magnetic perturbations induced in 
the planes, as for spinless impurities. In both cases, one 
can wonder whether the scattering saturates at low T or 
whether strong localization proceeds before this satura- 
tion can be detected. In other words, the MIT might re- 
sult from the conjunction of uncontrolled disorder and of 
the enhanced sensitivity of the underdoped regime to this 
disorder. This would explain why the MIT is strongly de- 
pendent on the cuprate family. 
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FIG. 1. The measured in plane resistivity p, scaled by the 
lattice constant c and the number n of Cu02 planes to give 
the 2D resistance R2D = pn/c, is plotted versus temperature 
for different irradiation conditions. The arrows indicate the 
sequence of measurements: increasing irradiation for straight 
arrows and long annealing at room T for wiggled arrows. 
a)T12201 : the data down to 1.7K has been taken in a differ- 
ent set up. b)YBC06.6 : the extrapolated low T variation of 
R2d{T) is shown as a full curve 

FIG. 2. The low T increase SR2D of the resistivity is plot- 
ted for R2D of about 4.8 fefi/D for the two compounds. 
The different fitting procedures used to estimate SR2D are 
described in the text: for YBCOe.Gjfits with a + bT^ or 
a + hT + cT^{a). For T12201, fit with a + bT^ between 70K 
and 150K (b) or between 150K and 300K (c). 

FIG. 3. a) Plot of Act versus InT for different defect con- 
centrations in Tl-2201. b) These curves are scaled by adjusting 
To in equ. (§). 

FIG. 4. The low T increase of R2D is normalized to the 
defect contribution to the residual resistivity Rd in YBCOe.e- 
Except for the onset of superconductivity, the data are found 
to collapse onto a unique curve as long as Rd < 4kQ/0. 
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Fig.2 - F. Rullier-Albenque et al 
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Fig.3 - F. Rullier-Albenque et al 




Fig.4 - F. Rullier-Albenque et al 



